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1 Introduction 

The influence of quenched, random disorder on phase transitions is of great 
importance in a large variety of fields Q, ranging from experiments with 
absorbed monolayers [Q in condensed matter physics to conceptual questions 
in non-perturbative quantum gravity 

For pure systems exhibiting a continuous phase transition, Harris Q de- 
rived the criterion that random disorder is a relevant perturbation when the 
critical exponent of the specific heat of the pure system is positive, a purc > 0. 
In this case one expects that the system falls into a new universality class 
with critical exponents governed by a "disordered" fixed point. For a purc < 
the behaviour of the pure system should persist, and a purc = is a special, 
marginal case. In two dimensions (2D) this scenario has been confirmed by 
various methods for many different systems ||. In three dimensions (3D) 
extensive computer simulation studies have concentrated mainly on the site- 
diluted Ising model ||, |(| . 

If a pure system with a first-order phase transition is subject to quenched 
disorder, the transition is softened and may even turn into a continuous one 
[0. This is always the case in 2D (§]; for numerical verifications see Refs. [9- 
13]. In higher dimensions, a tricritical point may appear at a finite concentra- 
tion of impurities juj, separating "non-softened" first-order and "softened" 
second-order regimes |lE| ] . Numerically such a scenario has recently been ob- 
served for the 3D siie-diluted 3-state Potts model |ll|. Since the first-order 
transition of the pure version of this model is very weak [ fl7|| , however, the 
characterization of the tricritical point remained inconclusive. 

In this report we give an overview on recent results obtained from ex- 
tensive Monte Carlo (MC) computer simulations of the 3D 2-state (Ising) 
[|l8| and 4-state Potts [19-21] models with &07id-dilution. The motivation to 
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study the 4-state Potts model derives from the fact that, in the pure case, 
this model is known to exhibit a fairly strong first-order transition, such that 
a disorder-induced softening to a second-order transition would give clear 
support for the theoretical picture sketched above. Modeling the disorder 
by bond-dilution enables in the Ising case a test of the expected universality 
with respect to the type of disorder. Furthermore, for both models this choice 
facilitates a quantitative comparison with recent high-temperature series ex- 
pansions p3, E3| for general random-bond g-state Potts models. 



2 Model and Simulation Setup 

The 3D bond-diluted g-state Potts model is defined by the Hamiltonian 

-pH = Y,Kij$*i,^ ai = l,...,q, (1) 

(ij) 

where the sum extends over all pairs of neighbouring sites on a cubic lattice 
of size L 3 with periodic boundary conditions, and the couplings Kij are 
distributed according to the distribution 

p(K l3 ) = P 5(K l3 - K) + (1 -p) 5(K l3 ), (2) 

where K = J/ksT is the inverse temperature in natural units. The dilu- 
tion parameter p is thus the concentration of magnetic bonds in the system, 
i.e., p — 1 corresponds to the pure case. Below the percolation threshold 
p c = 0.248 812 6(5) [M one does not expect any finite-temperature phase 
transition since without any percolating cluster in the system long-range or- 
der is impossible. 

The model (Q), (f| was studied by means of large-scale MC simulations us- 
ing the Swendsen-Wang (SW) cluster algorithm |25|| in the regime of second- 



order transitions, and multicanonical simulations C6, 27 in the regime where 
the first-order transition of the pure 4-state Potts model persists, i.e., at weak 
dilution close to p = 1. Thermodynamic quantities were averaged over a large 
number of quenched disorder realisations, ranging between 2 000 and 5 000. 
The stability of the disorder averages has been checked by monitoring run- 
ning averages as a function of the number of random samples. In fact, some 
care is necessary because too small a number of disorder realisations would 
lead to typical values rather than average ones p8| , and these two values are 
different if the probability distribution over the disorder realisations exhibits 
a long tail. 



3 Results 

3.1 3D Bond-Diluted Ising Model 

The phase diagram of the 3D bond-diluted Ising model has been obtained 
numerically from the locations of the maxima of a diverging quantity such as 
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Fig. 1. The average magnetic susceptibility \l of the 3D bond-diluted Ising model 
versus K — J/ksT for several concentrations p and L — 8, 10, 12, 14, 16, 18, and 
20. For each value of p and each lattice size L, the curves are obtained by standard 
histogram reweighting of the simulation data at one value of K. 
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Fig. 2. The energy autocorrelation time r e versus the size L of the 3D bond- 
diluted Ising model with SW cluster-update dynamics for various concentrations of 
magnetic bonds p on a log- log scale. The pure case corresponds to p = 1. 



the magnetic susceptibility \l depicted in Fig. |[ We focused on \l because 
the stability of the disordered fixed point implies a negative specific-heat 
exponent in a random system p9| . The error in this quantity is hence typically 
larger than that in the susceptibility. 

For an accurate determination of the maxima of the susceptibility we used 
the histogram reweighting technique with A^mcs = 2 500 MC sweeps (MCS) 
and between N s = 2 500 and 5 000 disorder realisations. The choice of Nmcs 
is justified by the increasing behaviour of the energy autocorrelation time 
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Table 1. Evolution of the 3D Ising model susceptibility for p — 0.7 and L — 96 
with the number of MC sweeps, JVmcs, for different disorder realisations, Xj: an< A 
the average value over 2 500 realisations, \. 
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T e as a function of p and L. At the critical point of a second-order phase 
transition one expects a finite-size scaling (FSS) behaviour r e oc L z : where 
z is the dynamical critical exponent. From fits to the data shown in Fig. || 
we obtained values of z « 0.27,0.38,0.41, and 0.59 for p = 0.4,0.55,0.7, 
and 1.0 (the pure case), respectively. The critical slowing down thus weakens 
for the disordered model and becomes less and less pronounced when the 
concentration of magnetic bonds p decreases. The largest autocorrelation time 
observed for the disordered model was around r e « 9 for p = 0.7 and L = 96. 
For each dilution and each size we thus collected at least 250 effectively 
uncorrelated measurements of the physical quantities, i.e., -/Vmcs > 250 r e . 
On the other hand, for small p « p c it is necessary to increase the number of 
disorder realisations because of the vicinity of the percolation threshold. 

In order to check the quality of the averaging techniques, we first studied 
the stability of the susceptibility versus the number of MC sweeps involved 
in the thermal average. For the largest size considered our results are given 
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Fig. 3. Distribution of the susceptibility for the different disorder realisations of 
the 3D Ising model with a concentration of magnetic bonds p = 0.7 and L = 96. 
The (running) average value over the samples X is shown by the solid line. 
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Fig. 4. Phase diagram of the 3D bond-diluted Ising model compared with the 
mean- field and effective-medium approximations (^) and (^), respectively. 



in Table [l] for different samples as well as for the disorder average. With 
2 500 MCS, the accuracy of the results for a given sample is not perfect, 
of course, but the precision of the average over disorder is quite good on 
the other hand. The disorder average procedure has been investigated by 
computing the susceptibility Xj f° r different samples, 1 < j < N s . As can 
be seen in Fig. |[ the dispersion of the values of x is n °t very large because 
the fluctuations in the (running) average value disappear already after a few 
hundreds of realisations. 

The phase diagram as obtained from the susceptibility maxima for the 
largest lattice size is shown in Fig. H. We do not include the results from 
high-temperature series expansions fl22|| in this figure since they would just 
fall on top of the MC data. For comparison we have drawn, however, a simple 
mean-field (MF) estimate of the transition point, 

K^{p) = K c {l)/p, (3) 

with K c (l) = 0.443 308 8(6) (fo). This rather crude approximation only holds 
in the low-dilution regime, p > 0.8. On the other hand, the single-bond 
effective- medium (EM) approximation of Turban |3l|| , 



Kf M (p) = lu 



(l-Pc)e K °M-(l-p) 



P-Pc 



(4) 



gives very good agreement with the simulated transition line over the full di- 
lution range. Since K c (l) andp c are input parameters this relation is trivially 
exact in the vicinity of both the pure system and the percolation threshold. 

Due to the competition between different fixed points, one of the main 
problems encountered in previous studies of the disordered Ising model was 
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Fig. 5. Effective exponents (l/u) c s as a function of 1/L m i n for p = 
0.95,0.9,0.8,0.75,0.65,0.6,0.45, and 0.36. The error bars show the standard de- 
viations of the power-law fits. The arrows indicate the values oil/u for the pure 
B3] and site-diluted [Ba] 3D Ising models, respectively. 



the question whether one measures effective or asymptotic exponents. Al- 
though the change of universality class should happen theoretically for an 
arbitrarily low disorder, it can be very difficult to measure the new critical 
exponents because the asymptotic behaviour cannot always be reached prac- 
tically. Another difficulty comes from the vicinity of the ratios 7/1/ and Pjv in 
the pure and disordered universality classes. Indeed, for the 3D Ising model 
these values are: 



7/1/ = 1.966(6), P/v = 0.517(3), v = 0.6304(13) pure case |32j, (5) 



j/u = 1.963(5), P/v = 0.519(3), v = 0.6837(53) disordered case §§. (6) 

Recent field theoretical determinations of critical exponents in the disordered 
case are presented in Refs. [35), and for an excellent review of various ex- 
perimental, theoretical and numerical estimates in the last two decades, sec 
Ref. ||. Thus, from standard FSS techniques, the critical exponent v only 
will allow us to discriminate between the two fixed points. This exponent can 
be evaluated from the FSS behaviour of the derivative of the magnetisation 
w.r.t. temperature which is expected to behave as dlnfh/dK oc L 1 /". From 
this power-law behaviour, we have extracted the effective size-dependent ex- 
ponent (l/v)cS which is plotted in Fig. || against l/L m i n for different bond 
concentrations p, with L m ; n denoting the smallest lattice size used in the fits. 
We clearly see that in the regime of low dilution (p close to 1), the system 
is influenced by the pure fixed point. On the other hand, when the bond 
concentration is small, the vicinity of the percolation fixed point induces a 
decrease of \ jv below its expected disordered value. This is plausible since 
the percolation fixed point is characterized by 1/u s» 1.12 p3j. 
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3.2 3D Bond-Diluted 4-State Potts Model 

Let us now turn to the 4-state Potts model which exhibits in the pure case 
a rather strong first-order phase transition. In order to map out the phase 
diagram of the diluted model we considered all concentrations p in the interval 
[0.28, 1] in steps of 0.04 and determined again the locations of the maxima 
of the susceptibility for a given lattice size L. The resulting phase diagram is 
depicted in Fig. ^, where we show for comparison again the simple mean-field 
prediction d3j) and the effective- medium approximation (||), using K c (l) = 
0.62863(2) [jigfl . On the scale of Fig. [| estimates from high-temperature series 
expansions up to order 18 are hardly distinguishable, for a comparison see 



In a second step, the order of the phase transitions was investigated. Here, 
in order to satisfy our criterion Amcs > 250 r e , the number of MC sweeps 
had to be increased up to 15 000 — 30 000, which is rather large compared 
to the values used in the Ising case. A first indication is given by the FSS 
behaviour of the autocorrelation time r e at the transition point. A glance on 
the log-log plot of Fig. ^ shows a crossover around p = 0.80 from a power-law 
behaviour for strong disorder (small p) to a clear exponential behaviour for 
weak disorder (p ~ 1), as is typical for a first-order phase transition. In fact, in 
the latter case one expects r e oc exp(2cr 0( j£ 2 ), where the (reduced) interface 
tension a d parameterizes the free-energy barrier separating the coexisting 
ordered and disordered phases. 

In the first-order regime we performed multicanonical simulations and 
estimated the interface tension from 
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Fig. 6. Phase diagram of the 3D bond-diluted 4-state Potts model as obtained from 
MC simulations as well as the mean-field (dashed line) and single-bond effective- 
medium (solid line) approximations (0) and (f3), respectively. 
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Fig. 8. Probability density of the energy of the 3D bond-diluted 4-state Potts 
model reweighted to equal peak height for p — 0.56 (top left) and p — 0.84 (top 
right). Interface tension versus inverse lattice size (bottom). 
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Fig. 9. FSS behaviour of the effective transition points for three different dilutions 
of the 3D 4-state Potts model as derived from the susceptibility maxima. 



where P max is the maximum of the probability density reweighted to the tem- 
perature where the two peaks are of equal height, and P min is the minimum in 
between, see Fig. ^| The linear extrapolations of a a( i in 1/L in the lower part 
of Fig. |8| imply non- vanishing interface tensions only for p = 0.84 and above. 
For p < 0.76, <j d seems to vanish in the infinite- volume limit, being indicative 
of the expected softening to a second-order phase transition. The tricritical 
point would thus be located around p = 0.76 — 0.84, in good agreement with 
the estimate of p — 0.80 derived from the analysis of autocorrelation times. 

To confirm the softening for p < 0.76 we have performed a detailed FSS 
study at p — 0.56 with lattice sizes ranging up to L — 96 and the number of 
realisations varying between 2 000 and 5 000 ||l9[ . The choice of p = 0.56 is 
motivated by our observation that in this range of dilutions the corrections 
to asymptotic FSS of the effective transition points are minimal, cf. Fig. |j| 

The log-log plot for Xmax in Fig- shows that for this quantity the 
corrections to asymptotic FSS seem to become quite small above L = 30, 
and fits of the form a x L 1 ^ 1 ' starting at L m in > 30 yield j/i> — 1.50(2). 
Using the data for L < 30 only, on the other hand, we obtained perfect fits 
assuming percolation exponents, j/v « 2.05 [^ij, cf. Fig. Similarly, the 
FSS of the quantity (din m/eLftT )x max oc L x l v gives for L m j n > 30 an estimate 
of the exponent 1/v = 1.33(3), consistent with the stability condition of 
the random fixed point (1/v < D/2 = 1.5) p9j ]. The same procedure was 
applied to the magnetization rh oc L~@/ u , but here the associated critical 
exponent turned out to be not yet stable. We therefore also considered the 
FSS behaviour of higher (thermal) moments of the magnetization, (fj. n ), which 
should scale with an exponent n(3/v. The results for the first moments exhibit, 
however, again much stronger corrections to scaling than we observed for \ 
or dlnm/ dK , leading to quite a conservative final estimate of (3/v = 0.65(5). 
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Fig. 10. FSS behaviour of the susceptibility, of dlnrh/dK , and of the magne- 
tization at T^max for the 3D bond-diluted 4-state Potts model at p — 0.56 (the 
quantities have been shifted in the vertical direction for the sake of clarity). The 
scaling behaviour for small lattice sizes below a crossover length scale is presumably 
governed by the percolation fixed point. 

We nevertheless note that our results do not fit satisfactorily the scaling law 
2/3 /i/ = d — 7/V. The reason could be the strong corrections to scaling at the 
random fixed point which are hard to cope with for medium-sized systems 




4 Conclusions 

By performing large-scale Monte Carlo simulations we have investigated the 
influence of bond dilution on the critical properties of the 3D Ising and 4- 
state Potts models. In the 3D Ising case the universality class of the disordered 
model is modified by disorder but its precise characterization turned out be 
difficult because of the competition between the different fixed points which 
induce crossover effects, even for relatively large lattice sizes. 

Applying similar techniques to the 3D 4-state Potts model we obtained 
clear evidence for softening to a continuous transition at strong disorder, 
with estimates for the critical exponents of v = 0.752(14), 7 = 1.13(4), and 
f3 = 0.49(5) at p = 0.56. The analysis of both the autocorrelation time and the 
interface tension leads to the conclusion of a tricritical point around p = 0.80. 
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